Probing the mechanical unzipping of DNA 
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A study of the micromechanical unzipping of DNA in the framework of the Peyrard-Bishop- 
Dauxois model is presented. We introduce a Monte Carlo technique that allows accurate determi- 
nation of the dependence of the unzipping forces on unzipping speed and temperature. Our findings 
agree quantitatively with experimental results for homogeneous DNA, and for A-phage DNA we 
reproduce the recently obtained experimental force-temperature phase diagram. Finally, we argue 
that there may be fundamental differences between in vivo and in vitro DNA unzipping. 
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Separation of double-stranded DNA (dsDNA) into 
single-stranded DNA (ssDNA) is fundamental to DNA 
replication and other important intra-cellular processes 
in living organisms. In equilibrium, DNA will denaturate 
when the free energy of the separated ssDNA is less than 
that of the dsDNA. Because of the larger entropy of the 
flexible single-strand, this can most easily be achieved by 
increasing the temperature of the sample until the DNA 
melts, somewhat above body temperature. In living or- 
ganisms, however, DNA separation is not only thermally 
driven, but also enzymes and other proteins may force 
the two strands apart. 

Recent advances in single-molecule force spectroscopy 
and dynamical force spectroscopy (DFS) has made pos- 
sible the systematic investigation of force-induced sepa- 
ration of dsDNA at room temperature where dsDNA is 
thermally stable in the absence of an applied force 
Although, these studies have significantly enhanced the 
understanding of the mechanical aspects of DNA repli- 
cation and transcription in vivo, it is also imperative for 
further developments of technologies, such as polymerase 
chain reaction and DNA chips, to understand the re- 
lation between thermal denaturation and force-induced 
separation. An initial step in this direction was taken by 
Danilowicz et al. |2| who published an experimentally 
determined phase diagram for the denaturation temper- 
ature as a function of the applied force. This study 
showed, as theoretically predicted that the force re- 
quired to unzip the DNA decreases with increasing tem- 
perature. However the applied theoretical framework Q 
does not capture all details in the entire temperature 
range. 

The mechanical unzipping of DNA has also been a sub- 
ject of several theoretical studies, which have often con- 
centrated on macroscopical aspects by investigating ther- 
modynamic equilibrium conditions |J| . Simulating realis- 
tic dynamics is unfeasible since the time scales reachable 
in molecular dynamics simulations are orders of magni- 
tude smaller than in experiments In this work we 
present a simple and efficient numerical Monte Carlo 
(MC) approach, to describe the unzipping process macro- 



scopically as well as semi-microscopically. In particular, 
we provide a theoretical underpinning for the experimen- 
tal force-temperature phase diagram recently published 
by Danilowicz et al. 

For this purpose, we use a simple one-dimensional 
model of DNA proposed by Peyrard-Bishop-Dauxois 
(PBD) 0. This model has been demonstrated 0, H @ 
to describe the thermally generated large amplitude local 
fluctuations quite accurately, an aspect which is gener- 
ally ignored in the thermodynamic models. We will show 
that this model successfully describes many of the micro- 
mechanical unzipping properties of DNA and provides 
significant insights into the physical phenomena govern- 
ing unzipping. 

In a typical unzipping experiment, one strand of the 
DNA molecule is tethered, at the terminating base, to a 
fixed surface while the other strand is connected through 
a polymeric linker to a force probe, such as a laser trap 
or an atomic force microscope (AFM) cantilever. For the 
sake of simplicity, the force probe is usually treated as a 
linear spring whose elastic properties can be determined 
through the system calibration. The force required to 
keep the molecule extended at a given distance is deter- 
mined by measuring the deflection of the force probe. 
By pulling the molecule apart at a constant speed, the 
force-extension curves of the system can be obtained. It 
should be mentioned, however, that the force determined 
by such techniques cannot be consider a characteristic 
quantity of the DNA, because it strongly depends on the 
elasticity of the experimental setup, as well as on the 
pulling speed. However, by performing the experiment at 
a wide range of pulling speeds and temperatures, useful 
information regarding the energy landscape of the DNA 
unzipping process can be accumulated. 

Similarly to this typical experimental setup, we use the 
model schematically represented in Fig. In this model 
there is only one degree of freedom describing the relative 
displacement of each base pair (bp) from its equilibrium 
position. The hydrogen bonds that hold the two strands 
together are approximated by a Morse potential (dashed 
lines), while the stacking interaction between successive 
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FIG. 1: A schematic representation of the PBD model in a 
dynamic force spectroscopy experiment. 
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FIG. 2: Force-extension curves of a 300 bp homogeneous AT 
sequence for three different values of the probe stiffness ko. 
The inset shows the height of the force barrier, F max , as a 
function of fco. All curves correspond to averaging over 10 3 
Monte Carlo simulations at T = 3007^ and vq = 0.1 V\ 



base-pairs is described by nonlinear springs (solid black 
line). The model's first base to the left is fixed, while the 
complementary base is attached to a linear spring (force 
probe) of stiffness fco which moves at a constant speed 
vq. The force probed by the spring is therefore given by 
Hook's law, Vp ull = fco(yo — V 1 )' wnere Hi an d Vo are the 
displacements (see Fig. ^| of the first base-pair and the 
spring's opposite end point, respectively. Specifically, the 
potential energy of the model is : 
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The PBD model parameters are those determined by 
Campa and Giansanti Q. To compute the force ap- 
plied during the unzipping process, we performed Monte 
Carlo simulations: The force probe is moved to the right 
(Uo ~ * J/o + Ay) by the distance Ay after which N Monte 
Carlo steps are performed in order to compute the aver- 
age displacement, (yi) of the first base pair. The force 
will then be given by F — fco(yo — (j/i))- The more MC 
steps used for sampling, the closer to equilibrium the 
system is before the next Ay move occurs. Thus, in this 
Monte Carlo framework, the pulling speed can be de- 
fined as = Ay/N. However, it should be emphasized 
that no direct comparison to the real time can be made 
as this depends on the detailed implementation of the 
Monte Carlo method. Here, the unit of pulling speed is 
V = 1CT 4 A/MC steps. 

We first investigate the initiation of the unzipping pro- 
cess. In Fig. |3 we present the force-extension curve of a 
homogeneous 300 bp AT sequence, for different values of 
the stiffness, ko, of the force probe. A significant force 
barrier at 2 A extensions is observed for stiff force probes. 
This barrier gradually decreases with decreasing stiffness, 
and finally vanishes for fco ~ 16 pN/nm (see inset of Fig. 
2). A simple interpretation of this result can be under- 
stood by considering the energy landscape of the entire 



system (DNA and force probe) [l(j. Assuming, for sim- 
plicity, the unzipping of the first base pair only, for small 
fco, the total energy has two local minima separated by 
an activation barrier. The first minimum corresponds to 
the bound state of the base pair and the second to the 
unbound state. As the force probe moves, the energy 
landscape is increasingly tilted, and eventually the sec- 
ond minimum becomes the global minimum. Beyond this 
point the unbound state is the energetically most favor- 
able. In the simulation (and also in the experiment) only 
the transition between the bound and unbound state can 
be observed. States close to the activation barrier are dif- 
ficult to probe, since the probability decays exponentially 
with increasing barrier height. For large fco, however, the 
elastic energy of the force sensor becomes dominant and 
the total energy has only a single minimum such that the 
unzipping process is continuous with no inaccessible re- 
gions. In fact, for very stiff probes (> 1.6 x 10 5 pN/nm) 
the resulting force accurately represents the derivative of 
the potential energy with respect to the extension. In the 
PBD model, the force barrier originates from a combina- 
tion of the force needed to break the hydrogen bond and 
the force needed to overcome the entropic barrier of the 
stacking interaction 0,^3- However, it is important to 
notice that this barrier can be observed only for the un- 
zipping of the first base pairs. When the first base pair is 
unzipped the effective probe consisting of the actual force 
probe and the newly formed single strand becomes very 
soft and the observation of the force barrier is no longer 
possible. An experimental indication of the existence of 
this force barrier can be found in the work of Krautbauer 
et al. 13], where significant force barriers are observed 
at the initiation of the unzipping process. However these 
authors attributed this barrier to interactions with other 
molecules. The experimental ambiguity of this barrier 
may be attributed to two factors: First, the experimen- 
tal apparatus does not have a resolution of 2 A where 
the barrier exists. Second, the typical total stiffness of 
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the force probe and the polymeric linkage is very small. 
The question of whether this force barrier exists when 
enzymes bind the DNA molecules obviously depends on 
how stiff their interaction is. The development of more 
sophisticated experimental techniques able to accurately 
probe protein-DNA interaction, will shed light on to what 
extend the binding process is driven by enzymes or ther- 
mal fluctuations y| and supercoiling effects 0] • In what 
follows, we use ko — 1.6 pN/nm, corresponding to the 
value of the force probe stiffness between the lower limit 
of AFM cantilevers and the upper limit of a typical op- 
tical tweezer. 

One of the most common questions investigated by 
DFS experiments is the dependence of the unzipping 
force on pulling speed. In Fig. Eta) we present the force- 
extension curve of a homogeneous AT sequence for three 
different values of Vq, at T = 300^. It is seen that faster 
pulling leads to higher measured force, in agreement with 
experiments 0]. However, for small vq (~ 0.2V) the 
measured force remains practically constant during the 
simulation, indicating that the system is close to equi- 
librium. This force corresponds to the experimentally 
observed unzipping force. During the unzipping of long 
DNA molecules, a slight increase in the slope of the force- 
extension curve is also observed, corresponding to the 
elastic energy of the stretched single strand and force 
probe system 0. It should also be noted that for even 
slower pulling speeds the molecule unzipping can take 
place more easily. This occurs because there is sufficient 
time for the system to be stochastical ly d riven over the 
activation barrier at even lower forces |l5j ■ 

According to our numerical results, presented in Fig. 
Efb), the unzipping force of a homogeneous AT and GC 
sequence is 20 and 36 pN, respectively. These values are 
rou ghly twice what one would expect from experiments 
fl6lll7l ] . To understand the source of this difference one 
should recall that the PBD model was originally devel- 
oped in order to study unforced thermal denaturation. 
In this context, the choice of only one degree of free- 
dom, which describes the relative motion of a base pair, 
was sufficient. However, this picture does not accurately 
describe the single strand dynamics in the kind of exper- 
iments we are considering here. During unzipping, two 
dynamically uncorrelated single strands form, and so at 
the macroscopic scale of the experiment there is a small 
possibility for two long single strands to be completely 
recombined. In the PBD model, the single strands are 
always dynamically correlated. As a result, in our numer- 
ical simulations, the single strands lack entropy, resulting 
in artificially high unzipping forces. In reality, j unzipped 
base pairs produce two single strands each consisting of 
j bases. The two single strands should be considered 
as a series of springs with effective stiffness k/2j, rather 
than as in the PBD model where the effective stiffness 
is k/j. A simple way of imposing this effective entropy 
reduction on the PBD model is to reduce the strength 
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FIG. 3: Force-extension curves obtained through Monte 
Carlo simulations for (a) AT sequence at three different 
pulling speeds, (b) AT and GC sequences at vo = 0.2V and 
vq = O.IV^, respectively. Subfigures (c) and (d) are the same 
as (a) and (b) but using the improved description of the 
singe strand stacking interaction (see text for details). In 
(d) vq — V and vo = 0.16V for AT and GC sequences, re- 
spectively. Each curve is the average over 10 3 independent 
simulations. In all cases T — 300A', ko — 1.6 pN/nm and the 
number of unzipped base pairs is 300. 



of the stacking interaction, in the single-stranded region, 
to k/2. This requires the introduction of a displacement 
threshold yth beyond which this transition takes place. 
The choice of this threshold does not critically affect the 
result, since our interest lies in describing the problem at 
the macroscopic scale of the experiments. 

In Fig. 13 (c) and (d) we present the numerical re- 
sults obtained from this approach. The resulting rate 
dependence on the calculated force is exactly the same 
as in Fig. Eta), but the unzipping force of homogeneous 
AT and GC sequences is 9 pN and 17.5 pN, respectively, 
values that are in excellent agreement with the existing 
experimental results 0, ^3 • It is important to empha- 
size that in vivo unzipping is expected to require more 
force than the experiments indicate 01 ■ Existence of the 
force barrier and the fact that the unzipping experiments 
probe the locally required force through a long and flex- 
ible single strand linker yields a lower unzipping force. 
In nature, enzymes apply the force directly to the DNA 
lacking the flexible linker that significantly lower the mea- 
sured force |15|. Indeed, simulations of the unmodified 
PBD model indicate that an enzyme is required to apply 
a force of at least 20 — 36 pN to unzip the DNA molecule. 
It is clear that our model and numerical approach will be 
able to shed new light on fundamental processes, such as 
replication and transcription when applied to interpret 
single-molecule exprimental data, obtain for these process 

With the Monte Carlo technique and the single strand 
modification to the PBD model described earlier, we are 
now able to determine the force required to unzip DNA 
at a given temperature and compare directly with the 
experimental results of Danilowicz et al. |2J . These au- 
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FIG. 4: Force-temperature phase diagram the A-phage DNA. 
Squares correspond to the experimental results of Ref. 0. 
Circles represents our numerical results for fco = 1.6 pN/nm 
and vo = 0.16V . Each point is the mean value of the force- 
extension curve obtained by unzipping the first 300 base-pairs 
of the molecule for 10 3 different simulations. 



involved in the unzipping of DNA. With the described 
modification, we were able to quantitatively reproduce 
the experimental force-temperature phase diagram re- 
cently obtained for A-phage DNA. Evidently, the PBD 
model successfully encompass both the traditional ther- 
mal separation and the more recently investigated force- 
induced DNA unzipping, and therefore offers significant 
predictive power for in vivo situations as well as emerging 
technologies. 

We thank Dr. Steven Koch for informative discussions. 
Work at Los Alamos is performed under the auspices of 
the US Department of Energy under contract W-7405- 
ENG-36. 



thors presented the force-temperature phase diagram of 
the A-phage DNA at a constant force, compared with the 
prediction of a simple thermodynamic model. Figure Q] 
depicts our results obtained for the unzipping of the first 
300 bp of the A-DNA. Since the molecule under investi- 
gation is not homogeneous, the unzipping force at each 
temperature corresponds to the mean value of the force- 
extension curve. As shown, our results are in very good 
agreement in the temperature range of 20 — 50 °C We 
notice a difference between our prediction and the ex- 
perimental results for T < 20 °C. This low temperature 
behavior may be due to a configurational change 2} , or to 
time scales of experimental techniques. Our predictions 
are more consistent with the experimental results than 
those of the thermodynamic model. This improvement 
is most evident at high temperatures, with the thermody- 
namic model yielding significantly higher forces. This is 
mainly attributed to the fact that the spontaneous bub- 
ble formation for T > 40 °C 21], is not described by the 
thermodynamic model. The bubble formation, which is 
predicted by the PBD model ||, leads to a considerable 
reduction of the measured force at those higher temper- 
atures. 

In summary we have presented a Monte Carlo ap- 
proach to mechanical unzipping of DNA, within the 
framework of the simple PBD model. Our method al- 
lows us to investigate the dependence on unzipping speed 
and temperature. We demostrated and analyzed the ex- 
istence of a force barrier at the initiation of forced unzip- 
ping which has generally been ignored in experimental 
setups. Further, we found that the single strand, as it 
extends between the double stranded molecule and the 
force probe, causes a decrease in the measured force. We 
showed that a modification of the PBD model is neces- 
sary to appropriately account for this effect and match 
experimental results. However, we find that the unmod- 
ified PBD model reliably captures the local dynamics 
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